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ABSTRACT 

The  reflection  coefficient  for  the  basic  mode  in  a  widening^ 
straight^  two-dimensional  waveguide  is  computed  for  small  wave 
numbers  by  using  the  perturbation  method  with  the  electrostatic  case 
as  the  unperturbed  case.   Essential  features  of  the  treatment  are: 

(i)   the  use  of  conformal  mappings, 
(ii)   the  consideration  of  physically  unreal  modes ;  and 
(iii)   the  inversion  of  an  infinite  matrix. 
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I.   The  Physical  Problem 
1.  Structure 

The  problem  which  we  are  considering  is  a  two-dimensional  one 
arising  from  a  three-dimensional  waveguide  structure  in  which  the 
perfectly  conducting  surfaces  extend  from  -co  to  oo  in  the  direction  of 
the  z-axis  in  a  Cartesian  coordinate  system.   The  intersection  of  these 
conducting  surfaces  with  the  x^y-plane  is  given  by  the  six  lines: 


X  >  0,  y  =  qrt 

X  >  0^  y  =   -qn 

X  =   0,  qjt<y<7t 

X  =   0^  -qjt  >  y  >  -T^ 

X  <  0,  y  =  K 

X  <  0,  y  =   -rt 


(0  <   q  <  1) 


The  segment  of  the  y-axis  between  -qr:   and  qn ,   which  we  shall 
refer  to  as  the  aperture,  separates  the  waveguide  structure  into  two 
simple  regions,  I  and  II. 
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2.  Conditions  on  the  electric  field 


Let  the  electric  field  E  =  (E  .E  .E  ) ,  and  assume  that  the 

X  y  z 

time  dependency  is  given  by  e    where 


u)  =  kc,   k  =  wave  number^   c  =  velocity  of  light. 


Then  we  have  the  conditions 


2 
1)  E  ,E  and  E  satisfy  Au  +  k  u  =  0; 
X  y      z 

Se   Be   Se 

5)  the  tangential  component  of  S  vanishes  at  the  boundaries. 
We  shall  also  make  the  following  assumptions: 


l^-)  E  is  an  even  function  of  y 

y 


5)  E^  -  0. 


The  behavior  of  E  at  infinity  must  satisfy 


6)   ^^^    lE"  -  tE  e^^l  =  0  ,  E  -  0 
•'  X  -»■  -00  I  y     o     '      '   X 


7)  ^^"^    IE   -  E  e^^-  pE  e-^^l  =  0  ,   E  ->  0 
x->-+ooiy    o         o      '         X 

8)  //  (E   +  E  )  dx  dy  is  boimded  over  every  finite  region. 


We  shall  introduce  the  following  definitions: 

T  (condition  6)  )  Is  called  the  transmission  coefficient. 
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p  (condition  7)  )  is  called  the  reflection  coefficient 

E  e    is  called  the  incident  wave. 
o 

E  -  E  e^^  is  called  the  diffracted  wave, 
y   o 

We  shall  now  show  that  the  field  is  completely  determined  by 
E  alone.   Conditions  2)  and  5)  give 

y 

Se    Se 


This  shows  that  E  determines  E  up  to  an  arbitrary  function  of  x 
y  ^ 

only.   But  since  E  and  E  satisfy  the  reduced  wave  equation 
■^  X      y 

Au  +  k^u  =  0;,  the  same  must  be  true  for  both  sides  of  2').   Furthemore, 
from  6)  and  7)  we  know  that  E  must  approach  0  at  co.   This  fact  gives 
us  the  uniqueness  we  seek.   For  if  2'),  with  a  fixed  E  ,  gave  rise 
to  two  different  solutions  for  E  ;,  then  their  difference  could  only 
be  a  function  of  x  alone.   If  this  difference  is  D(x) ,  then 


^  D(x)  +  k^(x)   =  0 


dx^ 


and 


D(x)  -♦  0  at  CO, 

This  implies  D(x)  =  0. 

Thus,  in  what  follows,  we  shall  limit  our  attention  to  the  determination 
of  the  y-component  of  the  electric  field.   Our  first  task  is  to  write 
the  boundary  conditions  on  E  . 


From  l)  ve  have 


a)  ZXE  +  k^   =   0. 
7  y 


From  3)  we  have 
5E 


SE 


d)  E  (o,y)   =  0  qjt  <  |y|  <  rt 


Note  that  in  b)  and  c)  we  used  the  fact  that  E  =  0  everywhere  on  those 

BE^  "" 

"boundaries.   This  implied  that  -^ —  =  0;  there)  and  then  2')  was  used  to 


get  b )  and  c ) . 

Continuity  conditions  in  the  aperture  are 


e)  E^^(0^,y)  =  E^(0",y)  |y|  <  art 


And  from  6)  and  7) 


V   lim   |_     „  ikxi     _ 
I)  E   -  tE  e      =   0 

X  ->■-<»'  y     o    ' 


,  ^   lam   |_    „  ikx    „  -ikxi 
h)         E  -  E  e    -  pE  e       =  0. 
X  -^  +00  I  y    o         o     ' 


II.   Method 

The  first  step  tcfu^ard  finding  the  reflection  coefficient  at 
the  interface  is  to  expand  the  field  in  each  of  the  regions  (l  and  II) 
in  Foiirier  series.   The  modes  are  determined  on  each  side  by  applying 
all  the  boundary  conditions  except  the  matching  condition  at  the 
interface.   The  solution  is  then  assumed  to  be  a  series  in  these  modes 
with  constant  coefficients. 

Applying  the  matching  conditions  at  the  interface  gives  rise  to 
two  infinite  sets  of  linear  equations,  one  arising  from  matching  the 
fields, and  the  other  from  matching  the  normal  derivatives.   Using  these 
two  sets  of  equations  we  can  eliminate  one  group  of  coefficients^  leaving 
an  infinite  matrix  equation  for  the  other  group  and  a  scalar  side 
condition  involving  the  reflection  coefficient.   Solving  this  matrix 
equation,  then,  will  allow  us  to  apply  the  scalar  side  condition)  and 
this  will  yield  the  value  of  the  reflection  coefficient. 

The  matrix  equation  which  arises  from  the  matching  conditions  is 
derived  in  full  and  is  valid  for  all  values  of  the  wave  number,  k  <  1. 
The  coefficients  of  the  expansion  of  the  field  in  region  II  are  the 
unknown  vector.   The  equation  is  solved  using  a  perturbation  method.   The 
electrostatic  case  (k  =  O)  serves  as  the  unperturbed  case.   We  shall 
see  that  the  ability  to  solve  the  electrostatic  case  for  arbitrary 
right-hand  sides  provides  us  with  enough  power  to  determine  the  higher 
order  terms  in  the  electromagnetic  case  (k  ^f  O) .   Using  this  fact,  and 

the  scalar  side  condition,  we  shall  determine  the  reflection  coefficient 

2 

up  to  terms  of  order  k  . 
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The  one  remaining  facet  of  the  method  is  the  proof  of 
solvability  of  the  electrostatic  case  for  arbitrary  right-hand  sides. 
Since  the  equation  in  this  case  is  Laplace's  equation^  we  may  use 
conformal  mapping  as  an  aid.   We  map  the  whole  waveguide  onto  the 
infinite  strip^  using  the  Schwartz-Christoffel  formula.   This  enables 
us  to  make  use  of  the  fact  that  each  unknown  is  a  Fourier  coefficient 
which  can  be  written  as  an  integral  (across  the  aperture)  of  the  field. 
The  field  is  written  as  the  derivative  of  a  harmonic  potential.   In 
this  way  the  matrix  equation  may  be  solved  in  the  electrostatic  case 
for  a  specific  right-hand  side.   However^  this  same  technique  can  be 
used  for  more  general  right-hand  sides.   In  order  to  accomplish  thls^ 
we  introduce  into  the  original  problem  "unrealistic"  modesj  i.e., 
modes  which  do  not  die  out  at  -H».   The  N   unrealistic  mode  will  give 
rise  to  a  right-hand  side  consisting  of  N+1  non-zero  entries.   Thus, 
solving  for  all  such  unrealistic  fields  will  yield  a  solution  matrix 
whose  product  with  the  original  matrix  is  a  triangular  matrix.   The 
entries  in  the  triangular  matrix  can  be  computed,  using  the  conformal 
mapping;  thus,  the  original  matrix  can  be  inverted.   Therefore,  we 
shall  have  solved  the  electrostatic  case  for  arbitrary  right-hand  sides. 

III.   Infinite  System  of  Linear  Equations 

In  all  that  follows  we  shall  use  the  notation  that  the 
restriction  of  E  to  region  I  is  denoted  by  E   and  similarly  for  E   . 

y  y  y 

When  referring  to  the  entire  waveguide  we  shall  write  simply  E  . 
Inasmuch  as  our  ultimate  goal  is  to  apply  a  perturbation  techaique  to 
this  problem,  we  shall  assume  that  k  <  1  so  that 
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n  -  k  q  >  0    for  all  integers  n  ^j:  0. 


By  separating  variables  in  a)  and  applying  conditions  b)  and  h)  in 
region  I  we  have 


\   (x,y)  =  E^  [_e    +  pe    ]  +  2_,  ""n^      ''°^  a' 

n=l 


k 

.  n 
1 —  x 


•  \A2  ,  2  2 

=  1  VE 


where  k  =  i  Vn  -k  q    n=1.2,... 
n 


In  region  II  we  have  the  additional  condition  d)  placed  on  the 

modes.   So  we  must  introduce  a  set  of  functions,  0  ,  satisfying  a),  c),  d) 

th 
and  g) .   In  addition,  let  the  N   function  satisfy 


e')   0^(o"^y)  =  cos  -  y  |y|  <  qrt 


in  the  aperture . 

The  reason  for  condition  e')  is  that  the  wave  coming  in  from  the  right 
is  a  superposition  of  cosines  in  the  aperture. 

It  \ri.ll  be  necessary  for  us  to  have  the  0  explicitly,  so  we 

th 
shall  derive  them  here.   The  N   mode  must  satisfy  the  conditions 

indicated  in  the  sketch  below. 


-3-^=0 

ay 


dy 


0 


^,  =  0 


N 


^N  =cos-q-y 


^,  =  0 


'^N 


Separating  variables,  we   see  that  0^^  can  be  expanded  in  a  series  of 
the  form 


r(x^y) 


(n)   iinx 
a    e     cos  ny 
n 


n=0 


-1  V  n  -  k. 


where  i  =  { 

n 


(region  II) 


n  4  0 


n  =  0 


At  X  =  0  we  have 


0jj(o,y) 


(N) 

a   cos  iry  =  <i 
n 


n=0 


cos  -y    y  <  qjt 

q 


q.rt  <  y  <  IT 


Consider  the  case  N  4  0^  ^^^   mtiltiply  both  sides  by  cos  my  and  integrate 


a           =     —      /  cos  —  y  dy     =      0 
o 

(N)             2       r  N        ^ 

a            =     —  cos  nry  cos  —  y  dy 

m               It    J  ^            q.         "^ 
o 


Thus 


S^^   e  cos  niy  W=l,2,. 

m=l 


where 


qit 


Similarly^   for  N=0,   we  get 


ikx  ^^  •^•^   ^ 

(5)  0Q(x;,y)      =      qe  +     ^       ^mo  ^     "^     ^°'^  °^ 

111=1 


The  field  in  region  II  may  now  be  written 


W  Ej^(x,y)   =  ^  d^  0^(x,y) 

n=0 


Applying  condition  e) ^  we  have 


Ejl  +  p)+^  c^cos^y  =  2_^    ^m  ^m^°^^^         |y|  <  q« 
n=l  m=0 
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■faut^  by  definition^  0  (o,y) 


m 
cos  —  y 

q. 


y|  <  q.rt 


Therefore^ 


(5) 


E  (1+p)  =   d 
o  o 


c   =  d       n=1.2., 
n      n         '    ' 


Differentiate  (l)  -with  respect  to  x 


ciEy  (x,y) 


=  IkE 


ikx     -ikx 
;    -  pe 


1 
+  — 

q. 


n 
1  —  X 

1    q.     °- 

c  K  e       cos  —  y 

n  n  q 


n=l 


S0„(x.y) 


m 


Let  \1/  (x.y)      , 
*m  '  dx 


Condition  f)  gives 


ikE  (l-  p)  +  - 


n 


c  k  cos  —  y  = 
n  n     q 


m  m 


n=:l 


m=0 


Solving  for  the  Fourier  coefficients  on  the  left-hand  side, 


qn       CO 


(6) 


i^o^^-P)   =     i    /     X!    'mV^°^^^    ^y 


o      m=0 


(7) 


'n       rtk  i 
n 


qn 


o       111=0 


<i^  t„(o,y)    cos  -  y  dy 
mm  q 


n=l,2,, 
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Combining   (5)   with   (7)^ 


(8) 


a  rtk  i 


n 


qrt 


o     ia=0 


d       t^(o,y)    cos  ^  y  dy 


n=l,2,. 


From  (2)  and  (3)  it  can  "be  seen  that 


ikx     V""         ^-^r"" 
(9)       toC^^y)   =  il^q.e    "*■  ^  /   ^ro  ^r  ^     '^^^  ^ 

r=l 


ilf  (x,y)   =  i   >   i      0       e     cos  ry 
^m  X  I   r  rm 

r=l 


m=1^2^. 


Therefore , 


(10) 


(i  ^l^  (o^y)   =  ikqd  +  i 
m  ^m  o 


m=0 


m 


m=0      r=l 


S   i  cos  ry 

rm  r 


Substituting  into  (8), 


q_rt 


n     rtk 


kqd  + 
^  o 


m 


m=0     r=l 


S   cos  ry  >  cos  —  y  dy 
r  rm       I      1 


n-1,2,... 


(11)  d 


m 


m=0     r=l 


k   rm  m 
n 


n=l,2,, 
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Equation  (ll)  is  an  infinite  system  of  linear  equations  in  the 

unknown  Foiu^ier  coefficients,  d  . 

n 

A  scalar  side  condition  is  obtained  by  substituting  (lO)  into 

(6). 

qrt 

/  00  oo 

if..  ^"^      ^"^ 


ikE   fl-  p)    =  —  <|  ikqd     +   i     >       d         >       S        i      cos   ry  >  dy 

o     ^  m=0  r=l 


(12)  =   ikqd     +     ~       >        d  >         -  S        I      sin  raq 

^  o  qn     /    I      ra     /     I      r     rm     r 

ni=0  r=l 


It  can  be  seen  from  equations  (5)  and  (5)  that  the  transmission 
coefficient;  -x ,   is  given  by 

(13)  T  -   q(l+p), 

which  gives  a  simple  relation  between  the 

reflection  and  transmission  coefficient.   Using  tlie  value  of  d 

^  o 

found  in  (5)^  we  have  from  (12) 


CO  CO 

(lif)                 >  d          >         -^     S        sin   qjrr     =     kE    (l-p  -q-pq) 

/    I  m     /    I      qjtr       rm            ^                    o          *^      ^     ^^' 

m=0  r=l 


This  scalar  condition  will  be  used  to  evaluate  p  after  (ll)  has 

been  solved  for  the  d  's. 

n 
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IV,   Perturbation  Method 

In  order  to  use  the  perturbation  method  in  the  solution  of  the 
equations  of  the  previous  section^  we  must  establish  their  explicit 
dependence  upon  k.   Having  done  that^  we  shall  solve  the  electrostatic 

problem  (k=0)  explicitly  and  use  it  as  the  unperturbed  case. 

i 
Equation  (ll)  depends,  upon  k  only  in  its  term^  -r—  . 

n 


r 
n 

\ 

"             r 

\/r2.  k^ 

r 

n 

Vn^-kV 

\  ~    \/T-^      ',/77?      ""V    /A     .2 


.-i^flM. 


/i2  .2  2 

n  1/     I   2    2 
V  r    n 


Equation  (ik)    depends  upon  k  only  in  its  term^  £ 


r 


o 


i^  =  -iVr  -k   =-irWl-—  =  -ir  +  2^  + 

'    r 


Equation  (ll)  can  be  written  nicely  in  matrix  notation  as  follows: 


11^ 


^  =  (^Q^^2.'^2'"'^ 


T  =  a  matrix  of  elements  T   =  —  S 

ixr       n     m 


n.  =/  0 


=   0 


n  =  0 


qjt 


S     =     matrix  of  elements   S       =  —      /   cos  ry  cos  — y  dy 

m   rt  J      ''  a"'  ^ 


(as  previously  defined) 


U  =  matrix  of  elements  U 


nr 


\  r    n  / 


=0  n  or  r  =  0 

Now  we  add  equation  (5)  to  equation  (ll)  and  use  the  notation  that  1 
is  the  vector 


1 


=  {l,  0,  0,  0,  ...  } 


We  assume  that  all  summations  go  from  0  to  +  oo. 


(15)      ^  =  -((T  +  k^)S 


■  (d)  +  E^(l  +  p)  1 


If  I  is  the  identity  matrix,  we  have 


(16) 


|l  +  TS  +  k^sj  (d)  =  E^(l+  p)  T 


*     i 

Let  T  = S   sin  qjtr 

rm    qn  rm 


c=0,l,2,. 


*      1    * 
rm    „  2  rm 


2r 


=  0 


r=l,2,.., 


r=0 
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Th.en  equation  {ik)   becomes 


(17)       >   d   >     (t"*"  +  k^*  )  =  kE^(l-  p  -q-pq) 
^  '  /    ,   m  /    ,     rm      rm      o' 

m=0     r=0 
For  shorthand^  let  us  define  vectors  t  and  u  whose  components  are; 

CO 

t   =   >   T 
m     /    I   rm 

r=0 


u   =   >   U* 
m     /    I   rm 

r=0 


Then  (17)  can  "be  rewritten  as 


(l8)       (t  +k^u*)-d  =  kE(l-p-q-pq) 


Equations  (lb)  and  (18)  fozm  the  system  which  we  shall  solve. 
First  we  shall  solve  (lo)  using  the  perturbation  method  to  invert 
the  matrix.   Then  we  shall  use  the  solution  vector ;,  d^  in  (i8)  to 
compute  the  reflection  coefficient,  p. 

The  method  will  proceed  as  follows.   We  shall  show  that  the 
electrostatic  system 


ll  +  TSJ  (d) 


can  be  solved  for  arbitraiy  right-hand  sides,  R,  in  such  a  way  that 

t  •  d  =  0.   When  this  is  known,  the  fiill  system  can  be  solved  up  to 

2 
terms  of  order  k  by  setting 
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d  =   e  +  k  f 


Equation  (l6)  becomes 


|l  +  TS  +  k^us]   (e  +  k^)   =  E  (l  +  p)l 


Thus 


/l  +  Ts]    (e)    +  k^  |l  +  TSJ      (?)    +  k^  jusj     (e)    =  E^(l+  p)l 

-> 

We  nov  can  find  a  unique  vector^  e^  such  that 

[l  +  Ts]  (e)   =  1 

(which  implies  that  E  (l+  p)e  solves  the  electrostatic  system)  and 
t  -e  =  0.   Now^  having  found  e^  we  compute  ] US  [  (e)  and  solve 

{l  +  Ts}  (?)   =   -  |usj  (e) 

in  such  a  way  that  t  •!  =  0. 

Having  found  ^(=e +  k  i )  we  substitute  this  solution  vect 


or 


into  (l8) 

kE^(i-p-q.-pq) 


(t%.^*,.(?..2?,  ,  ^^^-^ 


O 


But  t  -e  =  t  -f  =  0. 


Therefore , 


IT  - 


2->*  -.     KE^d-  p-q-pq) 

k  u  ■  e  =  — T^  /., s 

E^(l  +  p) 


(19)  p   =   l-(l+l^*-e 


1  +(q+  ku  -e) 


This  equation,  gives  the  value  of  the  reflection  coefficient  for  small 
values  of  wave  niimber. 

V.   Electrostatic  Case 


We  have  the  system  of  linear  equations 
[l  +  TS]  (d  )   =  t 


and  we  wish  to  show  that  it  caxi  he  solved  for  arbitrary  right-hand 
sides  in  such  a  way  that  the  solution  vector^  d,  satisfies  the 
scalar  condition 


t   •  d  =   0 


To  this  end  we  first  Invert  the  matrix  il  +   TSi  .   The  method  we  shall 
use  is  to  produce  a  matrix,  L ,   and  a  triangular  matrix.  A,  such  that 


{l  +  Ts]  [d} 


A 


Since  A  is  triangular,  its  inverse  can  easily  be  computed j  thus,  the 
full  solution  is 


[l  +  Tsj 


■1  =  DA-1 
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First  we  shall  construct  the  triangular  matrix^  A.   To  this  end, 
let  us  consider  the  field  in  region  I.   It  has  the  form 

00        n 
(20)  Ey^(x,y)   =  2^    ^n  ^  '^  cos  ^  y 

n=0 


in  the  electrostatic  case.   If  we  relajc  the  condition  at  +oo  and  allow 
the  field  to  become  exponentially  large,  new  modes  may  be  introduced. 
We  shall  refer  to  these  as  "unrealistic"  modes.   Let  the  N   unrealistic 
mode  be  denoted  by  E   (x,y) 


00  n 

/   N  x.  1/    >     \  '    N   "  q^      n 

(21)  E^  (x,y)   =  2_^      ^n  ^      cos  -  y 

n=-N 


In  region  II  we  have 


(22)  E^^^.y)   =  2_j  ^n     i^n^^^y) 

n=0 


where,  in  the  electrostatic  case. 


^Q(x,y)      =      q  +  p>       ^mo  ^"^  ^°^  ^ 
m=l 


(25) 


0^(x,y)      =      ^       S^  e        cos  my  n  ^  0 


mx 
b         f 

mn 

m=l 

a0 


o^  sr 1 

,    /        \  n  \  nix  , ., 

ijf„(x,yj    =  ^—  =        >       m  S        e        cos  my  all  n 

n  ox  /      I  Tnn 


m=l 
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Equations  (21)  and  (22)  provide  us  with  an  infinitude  of  new  fields. 
Solution  of  these  will  yield  the  matrices,  D  and  A,  which  we  seek. 

Matching  (21)  and  (22)  in  the  aperture  x=0,  |y|  <  qrt 


N      n 

c   cos  — y 

n      q-" 


^n     ^n^°^y^ 


n=-N 


n=0 


This   implies   that 


K 


=      d 


N 


n 


n  >  N 


(24)      c       +   c 

n  -n 


W 


0  <  n  <  N 


c        =      d 
o  o 


n  /^      n  , 

since  cos  —  y  =   cos      -  —      y 


anti  0^(o,y) 


cos  ^y  ,         \y\    <  qrt 


Let  us  define 


n 


=    < 


^     o 


n  >  0 


n  =   0 


1^ 

0 

n  >  N 

c^ 
-n 

0  <  n  <  N 

2      o 

n  =   0 
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Then 


(■25)  d    =  c   +  c  ~  ,  for  all  n. 

^  ^  n       n      n 


We  apply  the  derivative  condition 


m  c   e      cos  —  y 


dx  q  /    I      m  q_ 

m=-N 


^ —   =     /    ^^   K(^^y) 


n   n 
n=0 


(26)       >    —  (c   -  c  )  cos  -  y  =   >   d  \lf  (o^y)      lyl  <  qn 
m=l  n=0 


Multiply  by  cos  —  y  and  integrate 

qn 

CO 

n  /  -    +>      2    r   \  '  ^  N  ,  ,    ,      n   ^ 
-(c-c}=  —       >   d\^  (o.y)  cos  -  y  dy 
q  ^  n    n^      qn  J    /  ,   m  ^m^  '-^ '  q  "^   "^ 

o   m=0 


Multiply  this  by  —  and   add  equation   (25) 


qit 

CO 

(27)  2c"     =     d^  +     >       d^    —      /     t    (o,y)    cos  -Y  dy 

n  n        /    I      m     rtn    J      ^m'    ^'' '  0. 

m=0  o 
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However,  it  is  easily  seen  that  (27)  is  merely 
(28)  {l  +  TS]  (d^)   =   2c^- 


where 


N   ,  N 


^  =(C<.-) 


,  -*-      ,1  N   N         N   ^  ^    V 


If  for  each  N  we  can  compute  c^  and  d^ ,  we  will  have  found 
the  matrices  D  and  A.   The  solution  vectors,  d,^^,  will  comprise  Dj  and 
the  right-hand  sides,  c  ,  will  comprise  the  triangular  matrix,  A. 

Thus ,  we  must  do  the  following  for  each  N 
(l)   Compute  the  vector  c 

(ii)  Find  the  solution  vector  d^ 

Problem  (i)  may  he  solved  as  follows: 

We  map  the  waveguide  structure  onto  the  strip  |v|  <  rt  in  the  w^u+iv 
plane.   The  Schwartz-Chi'istoffel  formula  gives  us  the  mapping  implicitly 
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(29)  z     =      log 


yq  e^+1  +   qye^+l 


2q 


1-1  -1 


w  /'        2 
e     I  1-q 


V  q  e^+1     +    ye^+l 


dw 
dz 


2  w  , 
q  e  +1 


Let    ^ 


e      ^, 


a  =   e 


Substituting  into  (29) 


(30)  ^     =      a 


V  q^+a     +  \/l+?J  "^ 


V\/?7^  +   q^l^  J       [l  -  q^j 


l-£ 

q. 


Equation  (jO)  shows  that  a   can  be  expanded  in  a  power  series  in  ^. 


The  field  E  can  be  written  as  the  real  part  of  the  derivative  of -the 


N 
potential _,  e 


Nw 


T  4.  ^     d   Nw 
Let  (p,^  =  -r—  e 
■iN   dz 


so  that  E  =  Re  ^, 


i^N 


We  assert  that  ^  can  be  expanded  in  a  series  of  the  type 


Nz 


(1-N) 


,(N)  3T^,(N)   ^   q   %     ^^(N)  ^ 


^(N)/r 


1  =  1 


-  E  ^r^ 


(N).i 


i=-N 
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This  form  follows  from  (30)  and  from 


■X      ,,  Nw  dw     ^^  -N  /1+a 
i  =  ^^   d^  =  ^^   /-^ 

Y     Q    +0 


The   residue  theorem  asserts  that 


^0 


C-'       =      2^       J        ir      d^  1<    i<N 


This  reduces  to 


\/,2_  ,  ,/:n^  r  q 


q  +a  +  V/l+cr 


2i 


This  verifies  that  c^  .^=  0^   i  >  N. 

Problem  (ii)  may  be  solved  by  using  the  fact  that  in  the  aperture 
each  component  of  the  vector  d  is  a  coefficient  of  a  Fourier  series 
in  the  aperture. 

(52)      d^  =  -^   f  E,^  cos  2  y  dy  . 

-qjt 

But  E^,  =  Re  "S^,  =  Re  Ne^^  -^  . 
N      -*-N  dz 

Writing  the  full  complex  integral  (and  dropping  the  "Re"  for  now), 


2h  - 


qrix 


m  2rtqi 


„,     Nw     dw 

N  e  -r— 

dz 


-qiti 


m                m 
—  z         : 

e^     +  e      ^ 


dz 


2rtqi 


rti  +log(— ) 

e        <  e 


-«i  +  log(— ) 


2f(w)  -  t^^^^,  , 

q   ^        +  e      ^         I  dw 


where  F(w)    Is  the   Schwartz -Christ off el  mapping  given  in  equation   (29' 


Using  the   substitution  \  =  e 


w 


.N 


N 


(35)   d'  =    Re  ^^. 


(1-q)-    n    H-1  +  - 


vv?^  +  ^  ^^j 


A  2ni 


m 


dX 


\/  = 


q  X+1    +     yx+X 


(q-1)^     .    N-1-2 
+    (l-q^)  Mx  'I 


.\/ 


q^X+1  +  \A+1  j 


i  Vq  1+1  +    q\/x+l 


2ni 


d\ 


The  path  of  integration  is   around  a  circle   of  radius  —  in  the  positive 


sense ^   as  illustrated: 
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X-PLANE 


It   remains  now  to   show  that   solutions   of  the   system 


(l  +  Ts]     (d)      =     R  satisfy 


•X- 


t      •    d     =     0 


The  reason  for  this  is  that  the  components  d  of  d  are  the  Fourier 

m 


coefficients  of  the  derivative   of  a  potential  function  Re  ^ 
the   aperture. 


dv 

Sy 


m 


(51+) 


t    •    d 


2qrt 


qjt 


■q.rt 


'-^^^ 


Sv 


But  V  is  constant  on  both  boundaries^  so  that 


(o.qrt) 


(o^-qrt) 


Sv 


=  0  there . 


Therefore^  t  •  d  =  0 
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VI.  Convergence 

It  has  been  tacitly  assumed  in  the  previous  computations  that 
the  convergence  of  the  various  infinite  series  is  good  enough  to  justify 
the  operations  performed.   For  example,  the  final  formula  for  the 
reflection  coefficient  (equation  (l9))  involves  an  infinite  series 
represented  by  the  dot  product,  u  •  e.   We  shall  show  here  that  u  •  e 
does  converge  and,  in  fact,  that  it  converges  absolutely. 

From  equation  (17)  it  can  be  seen  that 


q.3t 

OO  00 

(55)  u   •    e  =     >       e  > sin   qrtr  -      /      cos   ry  cos  -  y  dy 

m=0  r^l  o 


_i 
qjt 


e 
2     /    ,      m 

m=0 


qrt 


cos  ry  cos  —  y  dy 


r=l 


,2         J  ''  q 


We  shall  first  show  that  the  factor  in  brackets  is  less  than  a  constant 

times  —  for  m  large, 
m 

It  can  be  seen  that 

/   cos  ry  cos  —  y  dy  =   (-1)   sin  mq 


■^o  ^  2  m 

r  -~ 

q. 


;n  the  special  case  of  q  =  —  the  above  formula  still  has  validity,  as 
can  be  seen  by  a  passeige  to  the  limit.   Thus,  we  shall  show  that 
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1 

r 


1      .  2 
-   sm  qjtr 


r=l 


2  m 


behaves  like  —  for  large  m. 
m 


We  break  this  sum  up  as  follows: 


1 


r=l 


2  m 

r       rH — 

q. 


1  .    2 

sm  rnq 


m. 
2q 


r=l 


2            m 
J    r       X 


1  .    2 

sm  rrtq 


"3ni 
2q 


1  1           .2 

-— -     sm  ntq 

2  m 

r       r 


m. 
2q 


+  1 


_jL_ 

2 
r       r- 


1         .    2 
sm  rjtq 


m 


r  = 


'Jm 
2q 


+  1 


where    [xj   =   integral  part   of  x. 
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h  < 

1 

m 

>: 

r=l 

1 

2 
r 

=  \^ 

<     - 

1 

m 
24 

2q 

.7< 

2ci 
m 

CO 

h 

> 

>; 

r=l 

r=l 

=     k 


2  m 


1 

^>~^^^            ' 

1 

1-^1 

< 

^^-^ 

h 

m 

^^'"^          1 

y 

2q 

r  = 

+  1 

4  m 


Finally ^    |s,| 


2q 


1  .    2 

sm  r:n:q 


r       1-  ■ 


m 


r  = 


in 

2q 


+  1 


Let  e 


m 


r  —   ,    so  that 


2i     <     ^r     <     2i 


For  each  r, 


1  .    2 

sm  rrtq 


m 
r  — 

q. 


1.2.  . 

, — ,    sm    (mjc  +   arte    ) 

hr 


.    2 
sm  qite 
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This  expression  is  clearly  bounded,  independent  of  m.   Call  this  bound 
B.  We  have 


<  B 

3iii 
_2q_ 

s 

>  ; 

r  = 

m 
_2q_ 

+  1 

2  <  ^^^\2 
r        /  m. 

2q 


m  =  k-,  — 
3  m 


Thus,  we  have  shown  that  the  original  sum  is  less  than  (k_,  +  k^  +  k^  +k,  )  — 

J.   2   5   4  m 


as  m  becomes  laree. 


In  order  to  complete  the  proof  of  convergence  of  u  •  e  (eq^uation 
(35))^  we  must  demonstrate  the  absolute  convergence  of 


(36) 


1 

—  e 
m   m 


m=l 


Each  component  of  the  vector  e  is  a  coefficient  of  the  Fourier  expansion 
of  the  field  in  the  aperture. 


qjt 


(37) 


m. 


1_ 
qir 


m 


E  cos  —  y  dy 

y    q 


-qrt 


Inasmuch  as  E  vanishes  for  qir  <  |y|  <  :n:,  equation  (37)  can  be  written 
as 


(58)      e   =  —  /   E  cos  —  y  dy 
m     qn    J  y  q.  ^      -^ 
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2rtq.i 


Jti 


■rti 


dw 
dz 


m  m 

—  z  -  —  z 

e'l       +  e      ^ 


dz 


We  can  convert  this  to  a  Line  integral  in  the  w-plane  using  the 
Schwartz-Christoffel  mapping  given  in  equation  (29). 


(59)  ^^    =     2^i   1 


yq  e^+1  +   q\/e^+l 

q. 


dw 


\/  2  w  _         i/v     7 
Vqe+1+    Ve    +1 


an 


m 


(1-a)^ 


q        C 


1-q 


-w       l^V  q  e  +1  +  /e  +1 


dw 


The  path  of  integration^  C,  is  illustrated  in  the  following  sketch: 


w-  plane 


I  I 

I     TTi+log  —5 
TTi  I  q  +7ri 

— r— 
I 


L 


TTI   I  .     ,  I 

-TTi  +  log     — 2 

q 


—  TTI 
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The  first  integral. 


(ho) 


2rtgi 


2ni 


\/?7^+  v^Xl) 


q. 


dw 


will  converge  if  the  path  is  altered  in  such  a  way  that  Re(w)  <  0) 
so  we  shall  change  the  path  into 

w-  plane 
+7ri 


-TTI 


Along  the  bottom  path  we  have 


w  =  -rti  -  s 


as  s  ->■  00 


Therefore,  the  absolute  value  of  the  integral  is 


(l-q)2  r 


1      ,.         2/   ^'q 


2s    (n/T??^  .  q\/r?^) 

y.  


2m 


ds 


But,  for  s  >  0, 
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\fr^  <  (v^::??^.  ^\/ZF^  <  (y 


<   VVl-q^e"""   +  Vl- 


Thls  means  that  the  absolute  value  is  less  than  or  equal  to 


1   ^1   2/   ^'q 


/,   sin   00    m 


2m' 


2m 


V/l-q^e-^  +\/^ 


ds 


1      2  f^-l)i 


m   00 


\^ 


^   2m     m 
.2m s 

l-q2  1     ^   e   1  ds 


2rtm 


The  same  analysis  goes  through  on  the  upper  path  of  integration, 
and  this  shows  the  absolute  convergence  of  the  first  half  of (56) 

For  the  second  integral  (59)  we  can  break  up  the  path  as 
follows : 


W-PLANE 


+  7ri 


•<r 


TTi  -^  log  -^ 


-TTI 


^.      q^ 
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Along  C    ,  \r  =   -k±  +   s   and  we  have 


2ni 


,,    ,  1  /         -m(-ni+s)     V  ^  V  y 


2.ai    I  i-ci^  J         ^    '°^72  ly?^  ^  ^\/^ 


-s 


1^ r   e- 


\   ^ 


ds 


^'"'^^  iog-%       rv7:?^.,vC 


2:rq.i^l-q     J  ^^2 


2m 


Along  C^^  w  =  rti  +   s. 


2m 


1  ■   ,     -'-                                      ^                    -                              \     n 

^     2                             ( \/  2      -s      y/        -s     )    ^ 

1                             /    ^       -m(ni+s)      V''  ^                 * 

/,      xm      ,/          ^                         ^^                                         \  2m 

27tqi( 

^^.^2^'""''^      "                           (V,^-.- .  ,Vl-e-J 

(42)   =^ /  e  -"^  ' ;:-  ds 

2jtai      _ 


G^^^  ^v/^^)  ^ 


q 

e  '  'r^  as 


,  ./,  2\     i  log-        l\/?:?^.^vw^ 

2nqi      1-q  q  ^ 


But,  a  comparison  of  (i^-l)  and  {k2)    (along  C^  and  C^)  shows  that  they 
cancel  each  other  out  exactly. 

We  are  left,  then,  with  integrals  along  C  and  C,  .   On  Cj^  we  have 


3U  - 


w=rtl   +    s  0<s<  log  —  2tii 


(1,3)   I r     e-"^^'^^^^^     — <i3 


\S    >J  /  , , \  2ni 


..,J.-,^y'"'^  -^^  i/.^-e- ..i/i-e- 


ly?:?^.^/: 


2  m 


-  ,m+l    r  q    Vy  1-q      V   1-q  y 
-1)       /         -^- ds 


.2   T   1   2 

1-q     "   1-q 


Note  that  in  the  Interval  both  numerator  and  denominator  are  of 
absolute  value  1.   We  introduce 


\  /  1-e       .      \  /e   -q     da    1 

cos  a=   \/— 2  .    srna=^--f-,  -=-^ 

^     1-q  "   1-q 


(W) 


"  1-q  '   1-q 


0  <   a,p  <  I 


-s 
e 

^e- 

;     2    ,/       -s 
-q      Vl-e 

s 
e 

\/.= 

-H/l-q^e^ 
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>5)       7  =  ^  -  P  .  S  =  -  I  "^^  P 


q. 


°<^<iU-\ 


Thus^  we  have  that  the  integral  along  C,  is 


lo  — 
2jt(ii      ^ 


In  the  same  way  the  integral  along  C   can  be  shown  to  "be 


1 
log  — 

M  ^       (-l)""    r  "^       -2nil7  , 

o 


Combining  {kb)    and  (1|7)  yields  the  single  integral 


log  — 

q 

(1^5)      lli^ /     sin  (2m7)  ds 


nq 


o 


We  now  wish  to  change  the  variable  of  Integration  from  s  to  7. 
Equation  (k^)    shows  that  7  is  monotonic ;  but  it  also  shows  that 
near  7=0  (s  =  log  —5)^  the  integral  will  become  improper.  Writing 
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and  letting  w  be  a  fixed  number  between  0  and  uj  ^  we  break  up  (kQ) 
as  follows 


(i^9) 


rtq 


-sin(2m7)2 cot  P  dj  +   '   sin  (Sm/)  —  dj 


w. 


For  any  fixed  value  of  uj  >  0^  cot  3  is  an  analy-tic  function  in  the 

interval  w  to  w  .   Thus^  we  may  integrate  by  parts  to  show  that  the 

first  integral  is  of  order  —  . 

m 

For  the  second  integral  we  must  go  back  to  equation  (i+3)-   We 
are  near  the  point  s  =  log  —  ,    so  we  set 


t  =   s  -  lo£ 


and  1  +  R  = 


1-e 


-s 


1-q 


+  i 


f   -s   2 

e  -q. 


1 
q. 


1  +  R  can  then  be  expanded  in  a  power  series  in  t,  near  t=0.   Tliis 
computation  shows  that 


(50)   1  +  R  =  1  -  s. A,^       +  higher  order  terms  in  t"""'^ 
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But 


2mij 


(l.H)^. 


Therefore^ 


(51) 


7  =  -i  log  (l  +  R)  =  -iR  +  higher  terms 


Combining  (5I)  with  (50)^ 


(52) 


7  =  13;l* 


5/2^ 


higher  tenns. 


Then  7  '  =  h   t   +  higher  terms.   This  mapping  can  "be  inverted  in 
some  neighborhood  of  the  origin  to  get 


.-I/2      1/3   u.  V,   ^     .    1/3 
t  '   =  c  7  '   +  higher  terms  in  7  ' 


This  Implies  that 


(55) 


22/3   ^.  ^ 
t  =  c  7    +  higher  terms. 


Thus 


(5^) 


ds  =  dt 


■   _  1      1  ■ 

d^7  5  +  f  (7^) 


d7 


where  f  is  an  analytic  function  near  7=0. 


We  are  now  in  a  position  to  return  to  the  second  integral  of 
(kS) .      From  the  foregoing  discussion  it  can  be  seen  that  w  should  be 
chosen  small  enough  to  be  within  the  radii  of  convergence  of  the 
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various  power  series  we  have  just  used.      That  being  the   case^  we  have 


(55) 


sin(2niy) 
.1/5 


dr/     +      /      sin   (21117)    f(7        )    ^7 


w. 


CO, 


The  second  integral  is  easily  seen  to  be  of  order  — ,  and  the  first  can 
be  evaluated  as  follows.   It  is  well  known  that 


(56) 


V 
X 


=  r  (1-  v)  cos 


nv  1 


2    1-v 

y 


Y  >  0,      0<v<2 


Now 


00     00 


w     o 


By  (56), 


sin  (2m7)  , 
1/5  ^-^ 


=  const 


.2/5 
m 


Further, 


(57) 


w. 


sin  (2m7)  -, 

1/5    '^ 
7 


cos  (2m7) 
an7^/5 


1   r  cos  (2m7)  , 
"1   '^ 


The  latter  integral  now  converges  absolutely.   This  shows  that 


39  - 


expression  (57)  is  of  order  —  . 

-*  -* 

Combining  all  these  results  shows  that  u  -e  converges 

absolutely.   For  we  have  shown  that  this  sum  is  at  worst  like 


m 


m' 


1 
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which  is  absolutely  convergent. 

VII. Concluding  Remarks 

We  have  restricted  our  investigations  to  a  first  order 
approximation  of  the  value  of  the  reflection  coefficient.   Higher 
order  approximations  are  available,  since  we  have  the  inverse  of  the 
matrix  determining  the  electrostatic  case. 

The  method  of  solution  used  in  this  paper  may  be  tedious, 
but  it  is  sufficiently  general  to  be  applicable  to  a  wide  range  of 
waveguide  problems.   For  example,  it  can  be  used  in  finding  the  field 
in  more  complex  waveguide  structures  such  as  those  illustrated  in 
the  following  sketches : 


B 


in 


IE 


IT 


1+0  - 


In  each  of  these  problems  there  are  three  distinct  regions, 
and  this  will  give  rise  to  two  sets  of  reflection  and  transmission 
coefficients.   (At  first  it  might  be  thouglit  that  regions  II  and  III 
in  problem  B  could  be  combined  into  a  single  region,  but  elementary- 
analysis  shows  that  a  field  of  this  form  is  over-determined  and  can 
yield  only  the  trivial  solution.) 


A  problem  treated  by  Magnus  and  Oberhettinger ' 


.s  essen- 
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tially  equivalent  to  our  problem  in  the  case  where  q  =  — .   However, 
the  method  used  there  is  based  on  algebraic  identities  rather  than 
on  conformal  mappings. 
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